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We show that the Ginzburg-Landau expansion of the grand potential for the Bogohubov-de Gennes 
Hamihonian is determined by the integrable nonlinear equations of the AKNS hierarchy, and that 
this provides the natural mathematical framework for a hidden nonlinear quantum mechanical su- 
persymmetry underlying the dynamics. 



I. INTRODUCTION: THE BOGOLIUBOV-DE GENNES HAMILTONIAN 

The Bogoliubov-de Gennes (BdG) Hamiltonian ; 1,] is 



where the (in general complex) function A(a;) plays the role of an order parameter, a condensate, or a gap function, 
depending on the physical context. The BdG Hamiltonian is relevant for the physics of: (i) crystalline condensates 
1i B in the chiral Gross-Neveu and Nambu-Jona Lasinio models 0, 0, 0, B @; (ii) fractional fermion number 
loTfllll: (iii) the Peierls effect [H, [11]; (iv) polaron crystals in conducting polymers [13, [H, [H, [13; (v) the pair 
potential in inhomogeneous quasi-lD superconductors [3, [iM [13, [HI ; (vi) the order parameter for superconductors 
in a ferromagnetic field (2^ ; (vii) incommensurate charge-density waves |23l . [2^ ; (viii) phase structure of ultracold 
atomic gases [1^ . In this paper we discuss the rich mathematical physics properties of the BdG Hamiltonian and in 
particular its close relation to nonlinear quantum mechanical supersymmetry [26l . [27| . 

An important object for the analysis of physical models based on the BdG Hamiltonian is the Gorkov resolvent 

R{x;E) = {x\j^\x) . (1.2) 

The resolvent ()1.2p is clearly a 2 x 2 matrix. It encodes the spectral properties of the BdG Hamiltonian. For example, 
the spectral function (density of states) characterizing the single-particle spectrum of fermions in the presence of the 
condensate A(a;) is 

p{E) = -liRTTD,rc [R{x; E + ie)] , (1.3) 

where the trace is a Dirac trace as well as a spatial trace. The local density of states is defined without the spatial 
trace: 

p{x]E) = -luYTrD[R{x;E + ie)] . (1.4) 

TT 

The thermodynamic grand potential density is defined as 

i^{x)^-\ r dEp{x-E)\n(l + e-^^''->''^) , (1.5) 

where (3 — ^ \s the inverse temperature, and p, is the chemical potential. The grand potential is given by the spatial 
integral 

, (1.6) 

and the thermodynamic properties [entropy, density, free energy, . . . ] of the system can be derived from [s^ . 
All information concerning the condensate A(a;) is contained in the local density of states p(x;E), and hence in the 
resolvent R{x;E). 
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For a one-dimcnsional system [as is being studied here] , the resolvent can be expressed in Wronski form, using two 
independent solutions of Hip = Eip. By elementary arguments 0, [2l|, [2^, [3^ , this construction implies that R{x;E) 
must satisfy a 2 x 2 matrix first-order differential equation, known as the Dikii-Eilenberger equation: 



|:i?(.;i^).3 



E -A{x) 
A*{x) -E 



, R{x; E) as 



(1.7) 



It must also satisfy reality and normalization conditions [expressed here for E not in the spectrum, with suitable 
extensions to other energy branches]: 

dcti?=-i , i?t = i? . (1.8) 
For example, if A = il/ is a constant, then we have the familiar solution 

^A f, ^i) (1-9) 



2VA/2 - E^ V ^ 

and it is straightforward to check that this resolvent satisfies the Dikii-Eilenberger equation p.7p . 

For non-trivial A(a;) it is not so easy to find the exact resolvent, so various expansions have been developed. One 
such expansion is based on a large E asymtotic expansion of the resolvent: 

n=0 

The energy integrals in (|1.5p and (|1.6p can now be done, leading to an expansion of the grand potential density and 
the grand potential 



■0(a;) = ^an(T,^i)gn{x) , (1.11) 

n=0 

oo « 

* = ^a„(r,Ai)G„ , Gn^ dxg^ix) . (1.12) 



This expansion generalizes a result of Perelomov [3J| for the grand potential of systems with Hamiltonians of 
Schrodinger form, rather than of Dirac form, and which is related in turn to the KdV hierarchy. On the other 
hand, we show in the next Section that the Dirac case is naturally related to the AKNS hierarchy. The T and ji 
dependent coefficients an{T,fi) are known functions 0, [1]. The gn{x) come from the diagonal elements of the 2x2 
matrices f„(a;) in the resolvent expansion (jl.lOp : 

f„(x) ^ f ^-r'^A . (1.13) 



Substituting the expansion p.lOp into the Dikii-Eilenberger equation p.7p . we obtain the simple recursion equations 

fn = -^/^i+Ag„ , (1.14) 



5; = z(^A*/n-i-A/,:_ij , (1.15) 

ffo = 1 , (1.16) 
/-I = . (1.17) 

The initial values follow from the A ~ M (constant) case above. It is simple to generate these terms recursively. The 
first few off-diagonal terms are: 



/-I = , 
/o = A , 



--(A"-2|ArA) 



f2 

h - ^(A"'-6|ApA') , 

/4 = (a('^) -8|A|2a"-2A2A*"-4|A'|2A-6A*A'2 + 6|A|4A 
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h = 



32 



(a(^) - 10|ApA"' - 10(AA*' + 2A*A')A" - lOAA'A*" - 10A*'A'2 + 30|A|''A' 



(1.18) 



The first few diagonal terms are: 



50 
.91 



1 


1, 



h = 



53 
54 
55 
56 



(AA' 



A'A* 



3|A|4 + 3|A' 



A"'A* - AA*'" + A'A*" - A"A*' + 6|Ap(A*'A - A'A*)) , 
L (^A('^) A* + A*('^) A - (|A'|2)" + 3|A"|2 - 10|A|2(A"A* + A*" A) - 5(A*2A'2 + A^A*'^) + 10|A|' 
L (^(|A|2)(H _ 5(|A'n" + 5|A"p - 10|Ap(A"A* + A*"A) - 5(A*2a'2 + A^A*'^) + 10|A|' 



16 
1 

32 
1 

32 



(1.19) 



Note that in principle there are integration constants appearing in the 5„, but these all vanish by virtue of the 
normalization condition (II. 



These recursively generated quantities /„ and ^„ have many remarkable properties, as is discussed in this paper. 
For example, the functional variation of Gn is related to fn-2 as: 



5G„ 



(n-1) 



(5A*(a;) 



fn-2{x) 



(1.20) 



The physical reason for this remarkable relation is the following. The variation of the grand potential (jl.6p can clearly 
be expressed as 



= 2^ari(r,Ai) 



(SA*(a; 



ri=0 



(5A*(x) 



(1.21) 



On the other hand, the grand potential (|1.6p can alternatively be written Q, including the thermodynamical Fermi 
factor in the energy trace, as 



dE 



1 



27r e^'(-E-A') + 

We choose to work in the Dirac gamma matrix representation 



- Tyd^. In [7" {E - H)] 



7=0-1 



1 

1 



7 



-1 

1 



7 = era 



1 

-1 



Then the BdG Hamiltonian H can be expressed as 



ax 



1(1 -7') A(x) + i(l + 7^) A*(a;) 



(1.22) 



(1.23) 



(1.24) 



4 



Therefore, the variation of the grand potential with respect to A* can be written as 

(n-1) 



= ^a„(r,M)^^/„-2(x) . (1.25) 

n— 2 

Comparing (|1.2ip and (|1.25p we find the non-trivial relation (|1.20p . 

II. AKNS HIERARCHY, FINITE-GAP SYSTEMS AND HIDDEN NONLINEAR SUPERSYMMETRY 

The resolvent expansion functions dnix) and fn(x) found in the previous section form the backbone of the (station- 
ary, defocussing) AKNS hierarchy [2^, [s^, HH, HH] , a hierarchy of integrable nonlinear differential equations. The idea 
behind the integrability of the AKNS hierarchy can be expressed in terms of a Lax pair [IBl, S and H, such that 

[S,H] = . (2.1) 

There exists an infinite number of such S commuting with H. This is a statement of the integrability of the AKNS 
system, and also of the BdG Hamiltonian. The Gorkov resolvent R{x; E) is the generating function for this infinite 
sequence of conserved quantities. To see this explicitly, define the matrix differential operator of order (iV + 1) by the 
following polynomial expansion: 

AT -1-1 

SN+i = iY,fN+i-k{x)a:iH'' , TV e No . (2.2) 

fc=0 

Here the ru are the same quantities that appear in the expansion p.lOp of the Gorkov resolvent. Then the recursion 
relations (|1.17p imply that 

[f^+i_,a3,i?] = 2f fN+i-k\_^f fNk\^ ^ (2.3) 



/Ar+l-fe / V /j 



iN-k 



This then implies that 

-2ir^Ux) 

The AKNS hierarchy is defined by considering finite linear combinations of the 5*^, with coefficients Ck'- 

N+l 

Sn+1 = ^ CAr+i_fe Sk ■ (2.5) 
fc=0 

Then the A^*'' equation of the AKNS hierarchy is [S'tv+i, H] = 0, which can equivalently be written as 

7V+1 

^CAr+i_fcA = . (2.6) 
k=0 

For example, the first few equations in the hierarchy are 

AKNSo : -^A'-hciA = , (2.7) 

AKNSi : -i(A"-2|A|2A)-ci^A' + c2A = , (2.8) 

AKNS2 : ;^(A"'-6|A|2A')-cii(A"-2|A|2A)-c2^A' + c3A==0 , (2.9) 
o 4 z 



AKNS, 
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^ (a('^) - 8|ApA" - 2A2A*" - 4|A'|2A - 6A*A'2 + 6|A|% 

-fci^(A"'-6|A|2A')-c2i(A"-2|A|2A)~c3^A' + c4A = , (2.10) 
o 4 z 
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The TV = 1 case, AKNSi, is the original AKNS equation, which can also be called the complex nonlinear Schrodinger 
equation. Another way to see the relation between the Gorkov resolvent and the AKNS hierarchy is to realize that 
the Dikii-Eilenberger equation is a statement of the zero-curvature representation of the AKNS hierarchy [2^, [s^l ■ 

From a physical standpoint, the AKNS equations arise naturally in a Ginzburg-Landau approach to the BdG 
system, in which one makes an approximate expansion of the grand potential (|1.6p to some finite order, A'^ + 1: 

JV+l 

^Ginzburg-Landau = Y,MT,f^)Gk . (2.11) 

Then the Ginzburg-Landau equation, ^^.'^^■j ~ 0, to this order is precisely an element of the AKNS hierarchy: 

N+i (k-^) 

5]afc(r,M)^-^A-2=0 (2.12) 

k=0 

with the coefficients expressed in terms of the real-valued functions ai^{T,^). Recall that /_2 = ./-i = 0. 

As already mentioned, the elements of the AKNS hierarchy have some remarkable properties, and because of the 
relation between the BdG resolvent and the AKNS hierarchy, these properties turn out to be very useful in solving 
physical problems based on the BdG Hamiltonian. 

For example, suppose A satisfies the {N + I)*'' order differential equation AKNStv- It automatically follows from 
the recursion relations (|1.17p that A satisfies all the higher equations AKNSjv+fe, for k ~ 1,2, . . ., with a suitable 
choice of the coefficients Ck ■ This means that there is an infinite number of conserved quantities for any member of 
the AKNS hierarchy. Physically, this has the important implication that it is possible to make a consistent ■polynomial 
ansatz for the Gorkov resolvent in order to solve the Dikii-Eilenberger equation p.7p . as was done in Thus, if we 
make the polynomial ansatz for i?(x; E): 



N 



n=0 



i?(x;i?) = ^/?„(i?) ( >'(^\ ^r;^f ) (2.13) 



then R{x; E) satisfies the Dikii-Eilenberger equation (|1.7p provided A{x) satisfies the AKNSAr_i equation. Therefore, 
given a solution of AKNS at- 1, the exact resolvent can be written in a simple polynomial form (|2.13p . with the 
coefficients Pn{E) determined from simple algebraic relations. This fact was used in Q to solve the gap equation 
in the NJL2 model, using TV = 2 in which case the condensate function A(x) had to satisfy the AKNSi equation, 
otherwise known as the (complex) nonlinear Schrodinger equation . In fact, this polynomial form of the resolvent of 
the BdG Hamiltonian generalizes an analogous property of Schrodinger Hamiltonians. For example, it was known 
already to Hermite that for the Lame potential, V{x) — TV(TV -I- l)j/sn^(a;), with TV integer, the resolvent can be 
written as a polynomial of degree TV in V{x) [Tsj . 

Next we connect this discussion to the known finite gap solutions of the AKNS hierarchy [sil, [s^. First, we discuss 
the spectral properties of the explicit finite gap solutions, and then explain their relation to nonlinear quantum 
mechanical supersymmetry. The remaining sections will discuss details of explicit cases where A is trigonometric 
(non-periodic case) or elliptic (periodic case). 



A. Finite gap solutions to the AKNS hierarchy 

When A satisfies an equation from the AKNS hierarchy, the spectrum of the system (|l.ip acquires a finite-gap 
nature [sil . [stI . [ssj . Suppose we have a complex periodic A that obeys the AKNS n equation [this picture can be 
extended for the case of a quasi-periodic A; see Section IV below]. The spectrum of (|1.1|) is composed then of TV 
valence bands, [i?2/c-i, £-2*;], k = 1,...,TV, and by two conduction bands, (— cx),i?o] and [i?2Ar+i, 00), separated by 
TV -I- 1 energy gaps, 

a{H)^{^(x^,Eo]U[Ei,E2]U...U[E2N-i,E2N]'J[E2N+i,^) ■ (2.14) 

The 2TV-I-2 energies Ek, k = 0, ...,2N+1, correspond to singlet band-edge states, which are described either by periodic 
or by anti-periodic wave functions. The states in the interior of the allowed bands are described by Bloch-Floquet 
quasi-periodic functions, and any corrresponding energy level is doubly degenerate, corresponding physically to left- 
and right-going waves. In the infinite period limit, A transforms into a non-periodic function, TV valence bands shrink 
and transform into TV bound states; the quasi-periodic states of the conduction bands reduce to scattering states with 
energies continuously varying in the intervals {—oo,Eq] and [i?2Ar+i, 00), where all energy levels, except edge values, 
are doubly degenerate. 
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B. Nonlinear quantum mechanical supersymmetry of AKNSjv 



A double degeneration of the spectrum, accompanied by a presence of a finite number of singlet states, is a typical 
property of A/" = 2 nonlinear quantum mechanical supersymmetry [26l . [39l . |40| . |4l| . In order to identify its structure for 
the finite-gap Bogoliubov-de Gennes system, we first describe briefiy different realizations of = 2 supersymmetry 
in Schrodingcr systems with a second order Hamiltonian. 

Quantum mechanical supersymmetry in its simplest form (4^ . [43l refiects coherently spectral properties of a 
composed system, given by superpartncr Schrodingcr Hamiltonians = —cP/dx'^ + V^{x), V^{x) = W'^{x)zLW' {x). 
The latter obey intertwining relations of a Darboux transformation [i^ , 

AiH- ^n+Ai , A\n+^n-A\ . (2.15) 

Here Ai = ^o-^-^ = ^ + W{x) is a first order differential operator, whose kernel is a (generically, formal) zero-value 
eigenstate iJo{x) of the second order differential operator . If ipo{x) is a nodeless function, the supcrpotcntial 
W{x) = —-^lntpo{x) and both superpartner potentials V^{x) are regular. Factorizations = A\Ai and W+ = 
AiA\ reflect a non-negative nature of the spectra of both subsystems, which in a generic case are almost isospectral 
by virtue of intertwining relations (|2.15p . In non-periodic case, when neither of the states ipo and l/ipo (the latter 
being a kernel of the operator A\) is physical, i.e. is not a normalizable function, two subsystems have exactly the 
same positive definite spectrum. If one of the functions ipo or l/ipo is normalizable, it is a singlet zero-energy state of 
the supersymmetric system 7i = diag(Ti.^ In periodic case both states ipQ and I /ipo are physical when they are 
periodic, and the lowest zero-energy level is doubly degenerate [47j . 

Intertwining and factorization relations are rewritten equivalently in the form of a linear (Lie) M = 2 superalgebra 

[Qa,W]=0 , {Qa,Qb}^25abn , (2.16) 

generated by the matrix Hamitlonian TL and supercharges 

Qi-^f °,t , Q2 = *^3Qi . (2.17) 



-A[ 

The integral of motion T ~ a^, with = 1, plays here the role of a Z2-grading operator, 

[F,W]=0 , {F,gj = . (2.18) 

Eq. (|2.15p is generalized for the case of a Crum-Darboux transformation [45| generated by an intertwining operator 
An of order n > 1, which annihilates n linear independent (generically, formal) eigenstates of 'H~ . The partner 
Hamiltonians are related by Ti"*" = — 2(lnyV„(?/'i, V'n))", where W„ is the Wronskian of the states of the 
kernel of A^- If W„(a:) is a nodeless function, the partner Hamiltonians Ti.^ have the same regularity properties. 
Intertwining relations are reformulated, again, in the form of zero commutators of the supercharges Qa, constructed 
in the form (|2.17p with Ai substituted by A„, with matrix Hamiltonian Ti = diag(7i+, 7i~). In correspondence with 
Burchnall-Chaundy theorem [H, [43], the linear superalgebraic relation (|2.16p becomes a nonlinear one, 

{Qa, Qb} ^ 26abPn{n) , (2.19) 

where P„(7i) is a polynomial of order n in the Hamiltonian. For such a non-linear supersymmetry, of order n, the 
operator F = 0-3 is again a Z2-grading operator. The spectra of the two superpartner Schrodinger Hamiltonians 
may differ in < I < n energy levels. 

Some non-extended systems with second order Hamiltonian are characterized by a hidden supersymmetry of a 
bosonized nature [3^, Ho, |4l|. In this case 7i+ = , and the intertwining relation reduces to the commutativity 
of H~ with a Hermitian differential operator which itself is identified as one of the supercharges, Qi = Z. Another 
supercharge is Q2 — iFQi, where a Z2-grading operator F is now a parity (reflection), or a twisted parity operator [40| . 
A bosonized supersymmetry of a nonlinear form appears, in particular, in finite-gap periodic and reflectionless non- 
periodic second order systems, in which a nontrivial operator of a Lax pair is identified as supercharge Qi annihilating 
all the singlet energy eigenstates. 

Corresponding systems extended by means of a Darboux or a Crum-Darboux transformation are described by 
a more rich, tri-supersymmetric structure (sol . [sil . [s^ . The latter admits three alternative choices for Z2-grading 
operator, F = (T3, 7?., or TZa^, where 7?, is a parity operator, and will be discussed in more detail in the next sections. 

For any type of supersymmetry described above, infinitesimal supersymmetric transformations generated by a 
supercharge Q can be presented in a form 



6X = [6aQ,X] , Sa = ^T 



(2.20) 
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where da is a composition of a Grassmann parameter ^ and of a Z2-grading operator T. 

For the BogoHubov-de Gennes system, Hamiltonian is a 2x2 matrix first order differential operator. When (in a 
generic case complex) the function A satisfies the AKNSat equation, the system displays a non-linear supersymmetry 
of order {N + 1). To identify its structure, we consider the 2x2 matrix operator Sn+i in (|2.5p as one of the 
supercharges, which annihilates the whole set of {2N + 2) band edge states, 

SN+i^Pk = , H^k^Ek^k , A: = 0,l...,2n+1 . (2.21) 

Using a matrix analogue of the Burchnall-Chaundy theorem, one finds that the square of 5jv+i produces a spectral 
polynomial of degree {2N + 2)^, 

2N+1 

S%+i = P2N+2{H) = \{{H-Et,) , (2.22) 

A:=0 

associated with the hyperelliptic curve of genus N . We will show that for any particular choice of the function A 
satisfying the AKNSjv equation, there exists an integral F, unitary equivalent to a nonlocal operator TZcr^i, with the 
properties 

[i/,F]=0 , {5jv+i,F} = , F2 = 1 . (2.23) 

Identifying integral F as a Z2-grading operator, we treat Sm+i as one of the supercharges, Qi = Sn+i, and define 
another supercharge in the usual way, Q2 = iFQi. As a result we find that the system is described by A/" = 2 nonlinear 
supersymmetry, 

[Qa,H]=Q , {Qa,Qb} = '25abP2N+2{H) . (2.24) 

The supersymmetric transformations arc given by Eq. (|2.20p . 

Notice that the described general structure of nonlinear supersymmetry of the finite-gap Bogoliubov-de Gennes 
system is similar to the hidden bosonized supersymmetry of a finite-gap Schrodinger system [40|, |4l| . In both cases 
nontrivial integrals annihilate all the singlet states, all the other physical states are organized in energy doublets, and 
the corresponding Z2-grading operators are nonlocal. 



III. REAL DEFOCUSING HIERARCHY NS+ AND TRI-SUPERSYMMETRY 



In this Section we consider the real case A = A* of the BdG system (|1.1[) . As it follows from a unitary transformation 
corresponding to a global chiral rotation with constant parameter Lp 

U{ip)H{l\)U-\ip)^H{K) , Uiif) = e^'P'^' , A = e'^A , (3.1) 

a real case of A can be understood modulo U{1). 

When A = A*, the Hamiltonian anticommutes with 17^ = cr2. As a consequence, the spectrum is symmetric with 
respect to the energy reflection E —E. Then the energies of singlet band-edge states (|2.14p obey a symmetry 

E2N+i-,=-E, , i = 0,..,N . (3.2) 

For the following discussion it is convenient to present Hamiltonian (jl.ip in another form by applying a unitary 
transformation generated by a matrix 



2^-1-* I -i 
This amounts to a change of Dirac representation from (|1.23p to 

7° -<^2 , 7^ = -^0-3 , 7^ = -o'l • (3-4) 
The BdG Hamiltonian (|1.24[) with real A takes therefore a form 

+ 
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This coincides with one of the matrix supercharges (I2.17P of a supersymmetric quantum mechanicanical system |42| 
described by a matrix Schrodinger Hamiltonian Ti. ~ H^, for which A = W{x) is identified as the superpotential, and 
supersymmetric partner potentials are given by 

y±(x) = A2±A' . (3.6) 

In accordance with (|1.23p [and also (|4.6p below], the operator ij^ underlying the energy reflection symmetry of the 
real BdG system corresponds to the grading operator T = for the associated superextended Schrodinger system. 

If A satisfies the AKNSat equation, the number of gaps modulo 2 in the spectrum of the periodic BdG system (or, 
of bound states in the non-periodic limit) depends on the properties of periodicity (normalizability) of the functions 

where u+ ~ (1, 0)', ~ (0, 1)*, and t means transposition. Namely, if iV = 2m — 1, rii e N, and both functions (|3.7p 
are periodic (or, one of these functions is normalizable) , they describe two (or, one in a non-periodic case) zero energy 
physical states of the system p.Sp with even number of gaps in a spectrum. For N = 2m, m £ N, states p.7p arc not 
periodic (normalizable in the non-periodic limit), and level E = does not belong to the spectrum characterized by 
odd number of gaps. 
From the viewpoint of the associated super-extended Schrodinger Hamiltonian 

n^Hl , (3.8) 

the cases = 2m — 1 and N = 2m correspond to finite-gap systems with exact and broken Af ^ 2 supersymmetries, 
respectively. 

When real A satisfies an AKNS equation, the equality (|3.8p . that connects the BdG and extended Schrodinger 
Hamiltonians, reflects an intimate relation between the AKNS and Korteweg- de Vries (KdV) hierarchies (see for 
example HH). This is just because in this case the Schrodinger potentials satisfy equations of the corresponding 
KdV stationary hierarchies. Therefore, each component Ti,^ of the diagonal Schrodinger Hamiltonian Ti has a non- 
trivial integral of motion in the form of a differential operator from the corresponding Lax pair (7i^, Z^) of the 
KdV hierarchy. The diagonal matrix operator Z ~ diag{Z^, Z^) commutes with the matrix Hamiltonian, [H, Z] ~ 0, 
and leads to existence of a more rich, tri-supersymmetric structure [sil. Is^. which will be discussed in more detail 
below. The systems with exact and broken J\f — 2 supersymmetry reveal the tri-supersymmetric structure in different 
ways, and these two cases will be considered separately. 



A. Exact Af — 2 supersymmetry and tri-supersymmetry in the light of the nS+ hierarchy 

Let us start with the case of unbroken J\f — 2 supersymmetry, and consider the BdG system p.Sp given by a real 
odd function A{x), 

A(-x) = -A(x) , (3.9) 

satisfying the (2?7i)'^ order equation AKNS2m-i, 'ni G N. Since Hamiltonian Hth anticommutes with both 0-3 and TZ, 
the operator TZas is a non-local integral of motion, [HR,TZa3] = 0. The stationary 2m — 1-th order AKNS hierarchy 
is described by the Lax pair {Hjgi, S2m), where 5*2771 is a non-trivial integral of motion, [HR,S2m] = 0. In the present 
case it anticommutes with the trivial integral TZa^, {5*2771, T^-Ca} = 0. We can treat TZas as the grading operator, 
and S2m as a non- linear supercharge for the BdG Hamiltonian in accordance with relations (|2.24p and (|2.23p . where 
A^ = 2m-1. 

Because of the presence of the nontrivial integral S2ni in the BdG system, the associated extended Schrodinger 
system with parity-even Hamiltonian Ti. is characterized by a more rich than Af = 2 supersymmetry structure. Indeed, 
relation (|3.8p implies that the matrix operator 52,77 is also the integral of motion for the associated Schrodinger system. 
Denoting Hs. = X and 527n = y, we have 

[H,x] = [n,y] = Q , [^-,3^1 = 0, (3.10) 

i.e. X and y are commuting parity-odd and parity-even integrals, respectively. As we noted above, the extended 
Schrodinger system Ti. is characterized also by the diagonal higher order matrix integral Z. Its relation with the 
original BdG system can easily be established. Consider the next after 52777 integral 52777+1 from the infinite set of 
integrals of the corresponding AKNS hierarchy. For the BdG Hamiltonian Hr, operator 5*2777+1 appears as a trivial 
composite integral. 



5277i+l — Hs,S2m — S2niHs, 



(3.11) 
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For the associated extended Schrodinger system, this is exactly the nontrivial diagonal parity-odd integral, S2m+i = 2, 

Z^Xy = yX . (3.12) 

Then relations ([3?T0l) can be supplied by 

[H,Z] = , [Z,X]^[Z,y]^0 . (3.13) 

Operators X, y and Z constitute a part of nontrivial integrals for the Schrodinger system. Hamiltonian Ti. commutes 
also with mutually commuting trivial integrals cts, TZ and TZcrs, and their products with X, y and Z produce other 
nontrivial integrals of motion. All these operators together generate the tri-supersymmetry of the second order finite- 
gap system Ti.. The resulting complete nonlinear superalgebraic structure depends on the identification of one of 
the trivial integrals, a^, TZ, or TZa^, as the grading operator. For any of the three possible choices, exactly one of 
the nontrivial integrals X, y, or Z is identified as Z2-even operator, which plays the role of the central charge of 
the superalgebra, while other two are treated as Z2-odd supercharges. This is in contrast with the supersymmetric 
structure of the BdG system iJu, for which the only imitary operator TZa^ can be chosen as the grading operator. 
For the details on the tri-supersymmetry we refer the reader to (sil . Is^ . 

Next we discuss some concrete examples of solutions and the associated explicit form of the nonlinear supersym- 
metric structure for the real AKNSat systems with N = 2m — 1. to S N. 



1. Real periodic solution of the (2m — l)th nS+ equation 

Consider a periodic solution of the (2to — 1)"^ nS+ equation (70j . 

^^^^ sn(..)cn(x,.) 

an [x, v) 

where to S N, and sn(a;, ^') = snx, cn (x, i^) = cnx and (li\.{x,v) = dncc are the Jacobi elliptic functions of elliptic 
parameter v, with < < 1 [ssj . The functions cna; and dux are even, while snx is odd, and therefore (|3.14p is an 
odd function. In this case the Hamiltonian (|3.5p takes the form of an odd operator 

d 





-j- TO Z -| 

ax an x 

Therefore, the system is characterized by a nonlocal integral of motion 



F = 7^cr3 • (3.16) 

Hamiltonian (|3.15p has a real period T ~ 2K(;^), Hs{x) — H]^{x + T), where K.{i') is the complete elliptic integral 
of the first kind. This system has 4to band edge-states of non-degenerated energies Ei, i = 0, 4to — 1, distributed 
symmetrically [see Figure [1]], 

E^ni-i-i = -Ei , i = 0, ..,2to , (3.17) 

with respect to the doubly degenerate zero energy level, which is in the middle of the central allowed band. All the 
band-edge singlet states are anti-periodic, their period is 4K(;^). 

With A satisfying the (2to — 1)'^ nS+ equation, the Hamiltonian (|3.15|) commutes with a non-trivial integral of 
motion (|2.5p . which can be presented in a form 

S2,n^( 1 , (3.18) 



where 



'^2m - „„2m4-l^ ( J„ „ J„ ) „„2m-l^ ' ^2m ~ i^27n)^ ■ (3.19) 



dn^+'x /cn^x^V" dn™-ix 
cn^^+ix \dnx dx ^ 

Since S'^^j is an even differential operator, p.lSp anti-commutes with the integral (|3.16p . {S2rmTZa3} = 0. Operator 
p.l6p is identified as a Z2-grading operator, and integrals 

Ql=S2m , Q2=^7^(73Ql (3.20) 
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FIG. 1: Spectrum of (|3.15|l for m = 2 and with i/ — 0.9. Red lines are the singlet band edge energies. The black line represents 
the zero energy level, with respect to which all the spectrum is symmetric. 



play the role of the supercharges for the system (|3.15|) . They generate the N = 2 polynomial superalgebra 

[Qa,HR]=0 , {Qa.Qb}^25abPiUHw) , (3.21) 

where, in correspondence with (|3.17p . 

PUHu)^I[lZ'\H^-E,)^U'X'iHl'Ef) . (3.22) 

The associated with p.lSp supcrcxtcndcd Schrodingcr system given by the second order Hamiltonian Ti = 
describes, in correspondence with Eq. (|3.17|) . the m-gap self-isospectral Lame system [U [mI [s^I with band-edge 
energies Eq = 0, E, = i^L.i+j = > 0, j = 1, . . . , 2m. 

For the BdG Hamiltonian in the simplest case m = 1, the singlet band-edge states and their energies are given by 

, / icn{x + b/2)\ ( isT\{x + b/2)\ ,,11 ^Q^ 

(,-cn(x-6/2)j ' ^^1= (,-sn(.x-V2)j ' ^^ = ^3^1 , ^3 - ^aV^o , (3.23) 

H^i,, = E,^, , Eo^-l , Ei^ -VT^ , E2 = -El , E3 = -Eo . (3.24) 
The explicit form of the states in the interior of the allowed bands can be found, for example, in refs. [s^l- 



2. Real non-periodic solution of the (2m — l)th nS+ equation 

Taking the infinite period limit of (j3.14[) . which corresponds to v 1, we get a non-periodic solution of the 
(2m — l)th nS+ equation, the familiar kink 

A^mtanhx , mGN . (3.25) 

Hamiltonian (|3.5p takes a form 

where the operator Vx satisfies the relations = — and V^^Dx = Vx+iD-x-i + 2A + 1. 

The central band [E2m-i, E2m], E2m-i = —E2m, of the periodic system shrinks here into a zero energy level, 
[E2m-i, E2m] = 0, and the corresponding band-edge states transform into a single bound state. The other 

valence bands shrink as well, producing (2m— 1) bound states of non-zero values of energy: [Ei,E2] £i, [E^, E4] 
£2, etc. The conduction bands transform into scattering bands with singlet edge energies £q < and £2™ = —£q > 0, 
where Eq £0, Ei^n-i = -Eo £2ra- 



11 



Defining the higher order differential operators 

2?!a=2?-a2?-a+i...2?-a+/3-i , P = 2,...,X , = 1, Vl^ = V^x , (3.27) 

we can write easily the eigcnfunctions of p.26|) . The 2m + 1 singlet states acquire a form 



HR4'tAx)=±£„^-l^f^/x) , Srn^l - ^"1^ " (™ " 0' • (3-29) 

Here / = m corresponds to two singlet edge states of the two symmetric parts of the continuous spectrum, while 
I ~ 0, ..,771 — 1 give (2m — 1) bound states. For I = 0, wave functions ^J^oi^) ^moi^) coincide, tp^Q{x) = 
V'm,o(2;) = V'm,o(a;), 

^"'«(^)= (sech™z) ■ ^3-30) 

This function describes a singlet bound state with energy £„i = in the middle of the spectrum. The scattering states 
of the bands (— oo,£'o] and [£2m,oo) have a form 

^± (X) - (-i + )^S^n;,V--l^e^^'^^\ 



Hmi^^i.^ix) = £m:k'ip^,.fix) , H^^ij^.^ix) = -Sm-k^^i-iix) , £rn:k V '"^ + , (3.32) 

where /c > 0. Boundary value fc = reproduces singlet edge states p.27p with I = m. 
The non-trivial integral p.lSp transforms into 



V_„,+iV_„,+2-V„,.iV„, 

nm \ p2m oil I?_„2?„,„+i...2?,„_22?,n-l 



(3.33) 



The supercharges are identified as in p.20p . They annihilate all the 2m + 1 singlet states. The nonlinear superalgebra 
is given by the same equation (|3.2ip as in the periodic case, but with the spectral polynomial operator changed 
coherently with the bands shrunk, 



m — l 2m — 1 



4=1 



The relation 5*1^ = PimiH^.) corresponds to a degenerate hyperelliptic curve. 

The associated second order supersymmetric Schrodinger system describes two reflectionless Poschl- Teller systems 
with m and (m — 1) bound states in their spectra. At m = 1 one of the two Schrodinger subsystems is a non-relativistic 
free particle [7l| . 



B. Broken M = 2 supersymmetry and tri-supersymmetry in the light of the nS+ hierarchy 

Now we consider the BdG system p.Sp defined by real A satisfying the AKNSat equation with N = 2m. As we 
noted above, in this case the associated extended Schrodinger system (|3.8p is characterized by the broken Af = 2 
supersymmetry. Our discussion will be restricted by three-gap BdG systems corresponding to the simplest case of 
m = 1, for which the explicit solution of the AKNS equation can be found in both periodic and non-periodic cases. 
This will provide us, particularly, with a new structure of tri-supersymmetry, to which the M ~ 2 supersymmetry of 
the associated one-gap Schrodinger system is extended. One of its peculiar features is that in periodic and non-periodic 
cases the super-partner potentials are given by the same function shifted for an arbitrary distance h. 

In the examples we discuss below, the nonlinear supersymmetric structure of the BdG system is based on the 
existence of the non-trivial local integral ^3 and a trivial nonlocal integral F = Tla2 playing the role of the Z2-grading 
operator. The first integral is a consequence of the AKNS2 equation satisfied by A, while the second means that it 
is an even function, A(— x) = A(x). Applying a chiral rotation, the latter integral can be reduced to the form of the 
grading operator we used in previous examples, i. e. TZa^. We do not apply, however, such a rotation here but work 
with the BdG Hamiltonian having the fixed off-diagonal form p.Sp . 
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From the viewpoint of the associated extended Schrodinger system, the non-trivial diagonal 2x2 matrix integral 

2 = 53 , (3.35) 

which is the third order differential operator, is related with the stationary KdV hierarchy given by the Lax pair 
[Ti, Z), \H, Z] = 0, where Ti. = X^, X = Hr. At the same time the next after S3 trivial integral 5*4 = S^Hr — HrSs 
for the BdG system is also a trivial integral for the Schrodinger system, 5*4 = ZX = XZ, cf. (|3.11[) and (|3.12p . 

Let us construct the matrix second order differential operator ^2 according to (|2.2p . A direct check shows that, as 
it should be expected, this operator is not an integral of motion for the BdG system, 

[Hm,S2] = ie{b)a3 , (3.36) 

where e{b) is some function of the displacement parameter, see Eqs. p.53[) and (|3.66p below. But since X = 
H]gi anticommutes with 173 (remember that this operator is a usual linear supercharge for the associated M = 2 
superextended Schrodinger system Ti., for which (73 is the grading operator), we get the relations 

[H, X] = [H, y] = [W, Z] = , (3.37) 
[Z,X]^[Z,y]=0 , (3.38) 
[X,y]^te{b)a3 , (3.39) 

where 3^ = ^2. Therefore, as in the case of the exact M ~ 2 supersymmetry, the three operators X, y and Z are 
nontrivial integrals of motion for the associated Schrodinger system. With respect to the trivial integral (73 identified 
as the Z2-grading operator F, integrals X and y are fermionic supercharges of the tri-supersymmetric structure, while 
Z is its bosonic central charge. Unlike the case of the exact Af = 2 supersymmetry discussed above, here the fermionic 
supercharges do not commute, see p.39p . Nevertheless, in both cases they satisfy a relation 

{X,y} = 2Z . (3.40) 

Since the BdG Hamiltonian iJR anticommutes with 173 and TZai, these operators alongside with TZa2 are the trivial 
integrals of motion for the associated Schrodinger system Ti. Therefore, besides the usual identification F = (73, 
nonlocal operators TZai and 7?,i72 can also be chosen as Z2-grading operators for the associated Schrodinger system. 
Under the choice F = TZcfi, the operators X and Z are identified as fermionic supercharges, while y plays the role 
of the bosonic integral for the system Ti.. For alternative choice F = TZa2, Z is again a fermionic supercharge, while 
the (boson-fermion) nature of the operators X and y is interchanged. For these two choices of the grading operator 
the commutator (|3.39p is a part of a complete list of (anti)-commutation relations of the tri-supersymmetry of the 
Schrodinger system. 

Let us stress that the essential difference of the described tri-supersymmetric structure of the Schrodinger system Ti. 
in comparison with the tri-supersymmetry from Section [ill Al ls that here the three trivial unitary involutive integrals 
of motion TZai , TZa2 and (73 anticommute mutually and satisfy between themselves the same set of algebraic relations 
as the three Pauli sigma matrices. At the abstract level, they generate non-abelian quaternion group of order eight in 
contrast with the Klein four-group generated by trivial integrals of the tri-supersymmetric structure of the extended 
Schrodinger systems associated with the family of the finite-gap BdG systems discussed above. 



Real periodic solution of the second nS+ equation 



Let us consider the periodic function 

. , , cn6dn6 , , , , , 

A(x) = \-usnbsn(x-b/2)sn{x + b/2) 

snb 

= Z{b-v)+'^^^^ + Z(:x-b/2-u)-Z{x + b/2-y) , (3.41) 
sn 

where Z{x]v) is the Jacobi Zeta bmction, and the real elliptic parameter is chosen, Q < v < \. The real parameter 
b lies in the range < 6 < K(i^), where K(i/) is the elliptic half-period. The function (|3.4ip is a solution of the 
third-order AKNS2 equation (to be more specific, in the terminology of [sij , the second nS-(- equation): 

A"'-6A2A'-2(l + z/-3ns2fe)A' = . (3.42) 
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The Hamiltonian (13.51) takes the form 



_d_ 

dx 



' cn b dn b 
^ sn 6 





+ i^sn6sn(a;-fe/2)sn {x + b/2)) 



_d_ 

dx 



+ 1/ sn 6 sn (.t - 6/2) sn (x + 6/2) 




(3.43) 

and we can check that it has a trivial integral of motion F = TZa2, [7?.(T2,Hr] = 0, where TZ means reflection. The 
square of (|3.43p produces a self-isospectral Schrodinger operator with potentials = ± A' as in (|3.6p , satisfying 
the translation property V~{x) = V~^{x — 6): 



^ + ns^ 6 - 1 - 1/ + 2i/sn2 (x + 6/2) 



dx 









+ ns^ 6 - 1 - z/ + 2:/sn2 {x - 6/2) 



(3.44) 



The spectrum of this system [see Figure [2] is composed of four bands, arranged symmetrically with respect to the 
zero energy level, which does not belong to the physical spectrum. 



Es 
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{a.)v = 0.2,6 = 1 



{c)u = 0.7,6 = 1 
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{h)u = 0.9,6 = 1 



(d)y = 0.7,6 = K(0.7) 



FIG. 2: Plots show the spectrum of p.43|) and its dependence on v and b. The red lines correspond to the band edge energies. 
Figures (a) and (b) show the cases with different values of the modular parameter v and the same value of b; in figures (c) 
and (d) the displacement constant 6 varies, but v is the same. The last plot (d) corresponds to the limit case when (|3.41[) is 
reduced to (|3.14p with m = 1. 



The six band-edge states can be written easily in terms of the Jacobi elliptic functions 
, / idn (x + 6/2) \ , f icn (x + 6/2) \ , 



_ f idn (x + 6/2) 
~ -dn (x - 6/2) 

■03 



^ icn {x + 6/2) \ 
^-cn (a; -6/2) ) 

, ■04 = 0-301 , = o'aV'o 



zsn {x + 6/2) 
-sn(x - 6/2) 



dnb 

sn b 



E, 



'E. 



cn b 
sn b 



(3.45) 
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From the real nS+ hierarchy we find a non-trivial integral of motion S3 which annihilates all the band edge states, 



S3^^{% g-] , [i?M, ^3] = , (3.46) 



3 

where 

d 

ii'dn {x + 6/2)cn {x + 6/2)sn {x + 6/2) , 
S^{x) = S+ix-b) . (3.47) 

The integral of motion 6*3 anticommutes with TZa2 , and we can identify the M = 2 nonlinear supersymmetry. It is 
generated by the supercharges 

Ql=53 , Q2 = i7ea253 , (3.48) 

[-ffE,Qa]=0 , {Qa,Q4 = 2<5,faP6(i?R) , a, 6 =1,2 , (3.49) 
where Pe{Hs) is a spectral polynomial 

5 

Pe (Hm) = n (^« - ^fc) = (^K - ^0 ) (^^R - ^^1 ) (^E -El) . (3.50) 
Another operator from the hierarchy which commutes with 6*3 but does not commute with i^R is 

S2=t(^ g+ ^ , [52,^3] =0 , (3.51) 

where 

^2^ = ^TA^-i(A^T A'-l-. + 3ns^5) . (3.52) 



We can check that 5*2 satisfies 



dn 6 cn 6 , , 

[HM,S2] = 2ia3 — . 3.53 



Thus, this example shows how the nonlinear supersymmetry structure arises for the finite gap BdG system given 
by the periodic function p.4ip . We note here that the Schrodinger system p.44p was also discussed in [sl] in the 
context of Darboux transformations. 

2. Real non-periodic solution of the second nS+ equation 

Now we consider the non-periodic limit of (|3.4ip when the elliptic parameter v 1, 

A{x) = coth6 + tanh(a;-6/2) -tanh(a; + 5/2) . (3.54) 

The main properties of the supersymmetric structure remain invariant in this limit. The gap function p.54p is also 
a solution of the second nS+ equation 

A'" - GA^A' + 2(1 + — ^)A' = . (3.55) 
sinh b 

The Hamiltonian (13.51) takes here a form 



_ Y ^ ^ + tanh(x - 6/2) - tanh(x + 6/2) , 

"'" ^i^ £-(coth6 + tanh(2;-6/2)-tanh(a; + 6/2)) ' ■^'^■^'^^ 
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The square of JTr produces two self-isospectral reflectionless Schrodinger Hamiltonians, 

The spectrum of p.56p is symmetric with respect to the zero level of energy. It is composed by two bound states 
which come from the shrink of the bands E2\ £i = ^ 5557175, and [i?3, £^4] £2 ~ 51375' 



(irecht + vl) ' ^2 = ^3^1 , H^i^., = E,^, , 1 = 1,2 , (3.58) 



and by two continuous bands (— oOj^o] and [£4,00). The energies in the continuous bands arc doubly degenerate, 
except the singlet edge values, and arc described by the eigenfunctions 



=Ffc+icothbp+ ±»fc3: ^ ^ ^ a/(1 + sinh^ 6 + 1 



, , ^ ^ 1^ \ + + y 1,1 + K-j smn^ / 

^t=[ ) ' Wf=^T.^f , SV,-^- , (3.59) 

= '^s^f , i^aV'f = ^i^fcV'f , ^Uk = -£v,k , (3.60) 

where we introduce the notation 



= + A tanh(a; ± b/2) . (3.61) 



The singlet edge energies £i-q = £4 and £i-o = £0 correspond to zero value of the non-negative continuous parameter 
k. 

The integral of motion 5*3 of the corresponding Lax pair (i?K, 'S'3), which annihilates all the singlet states, can be 
obtained by taking limit — > I in p.46p . It has the following structure: 

As in the periodic case we have 

[Hm,t] = o , {53,r} = o , r = a27e , (3.63) 

and the supersymmetry has the same non linear structure (|3.48p . (|3.49p with the spectral polynomial changed for 

Pe{HR) = {Hl-£^){Hl-£!f . (3.64) 
The components of the operator 5*2 from ()3.5ip reduce in the infinite period limit to 

St = -f^ T (coth6 + tanh(x - 6/2) - tanh(a; + 6/2))-f- H + sech^ {x ± b/2) . (3.65) 

dx-^ dx sinh h 

This operator obeys the relations 

[5'2, gg] = 2i '^.°^^3^ 0-3 , (3.66) 
smh 

[52,53] = . (3.67) 
IV. COMPLEX CASE 

Here we describe a one-parametric generalization of the finite-gap systems from Section IIII Al given by a complex 
A. In such a case an energy reflection symmetry is partially broken, but the supcrsymmetric structure of the BdG 
system is preserved in a form to be coherent with the partial discrete symmetry breaking. A kind of duality underlying 
the supcrsymmetric structure is observed. 
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A. Complex non-periodic solution of the (2m — I)*'' nS+ equation 

First, consider a simpler non-periodic limit case. We take the following complex generalization of the kink (|3.25p . 
the twisted kink [3, [s^l , 

^ ' ' cosh ( sin (6'/2) x) ' ^ ' 

where m G N, and 6* is a real parameter. The function A(.t) in (|4.1[) satisfies the complex {2m — 1)*'^ nS+ equation. 
It obeys relations 

A(x; 6* + 27r) = A(-x; 61) , A(a-, + 47r) = A(a;; 6*) , (4.2) 

and reduces to the real kink case ()3.25|) at 6* = tt. It takes constant values +m and — m at 6' = and 6 = 2-11 
respectively, when the BdG system reduces to a free (1+1)D Dirac particle of mass m. 

In correspondence with (|4.2p . the BdG Hamiltonian (|1.24|1 anticommutes with the following two mutually commuting 
operators 

^1 = -C/(f)7°T2,t/-i(f) , ^2 = C/(f)7^^2.C/-l(|) . (4.3) 

That is, {H, Qi\ = 0, i = 1, 2, where U is the operator of the chiral rotation defined in (|3.ip . and T^ti is an operator 
of the half-period translation of 9, T2-K : ^ 9 + 2tt. Notice that the treatment of the operators (|4.3p in complex case 
9 TT implies the inclusion into consideration of a "dual" BdG system with shifted in 27r value of the parameter 9. 
Because of the presence of the operator T2tt in (|4.3p . the reflection energy symmetry of the real case, as we shall see 
below, is partially broken here. On the other hand, the half-period translation T2tt disappears in the composition of 
operators (|4.3p . 

^ = 7^C/(|)7^C/-l(f) , (4.4) 

which is a nonlocal integral of motion of the system, [F, H] = 0. 

We present the BdG Hamitonian (|1.24p in a more suitable form applying a unitary transformation given by the 
matrix 



which reduces to (|3.3p when 9 = n. Equivalently, this corresponds to a choice of the following 6'-dependent represen- 
tation for 7-matrices, 

/ sin I cos I \ , / - cos | z sin | \ , 

^ =Ucosf zsinfj ' ^ =U-inf cos f j ' ^ ^ ' ^'''^ 

In this representation the BdG Hamiltonian (jl.24p acquires the form 

^ / imcos| ^ -f msin I tanh(sin |x) \ /— mcos| i2?,„ \ (A 7) 

\ ^ ~ "^^i'^f *'^^l^(^i^f^) -imcosf / V iP-™ mcosf/ 

In (|4.7p we introduced the same type of operators as in the previous section, 

d 9 f 9 \ 

i'x = — h A sin - tanh i sin -x j , (4.8) 

which satisfy the relation 

Vxi^x) = sin f 2?A(sin fx) (4.9) 
and reduce to I?a, defined in (|3.26p . when 9 = tt. The operators defined in (|4.3p take in representation (|4.6p the form 

Gi - <y3T2^ , Q2 = nT2^ , (4.10) 
and the nonlocal integral (|4.4p reduces to p.l6p . F = Tlcrz- 
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Hamiltonian (|4.7|) satisfies the relation T2-kHc = {He + 2m cos |(T3)T2t, that transforms into [T2^,iJc] = in the 
case 9 = TT. Its eigenvalues depend on 9 and satisfy the relation £{9 + Att) = £{9) in correspondence with (|4.2p . 
Since He anticommutes with operators (|4.10p . we find particularly that if a state is an eigenstate of (|4.7p of 

eigenvalue £{6), then a state TZT2-,r4'S{,e){x) is an eigenstate of eigenvalue —£{9 + 27r), 

Hcipei,e){x) = £{9)^£ie){x) => HcnT2^^e{e){x) ^ -£{9 + 2^)nT2^^£i^e){x) . (4.11) 

Explicit form of the eigenstatcs and eigenvalues of the BdG Hamiltonian (|4.7p . discussed in Appendix, shows that 
for 9 ^ Q (modTr), if £{9) is an energy of a physical state ^£(^g){x) to be different from the middle bound state, then 
the state crzT2TT'4^£{0){x) (or, TZT2Trip£i^g){x)) possesses the energy ~£{9). Together with (|4.1ip . this means that all the 
energy levels except the energy of the middle bound state are periodic functions of 9 with the period 27r. 

In more detail, as in the real case, for 7^ 0, 27r the spectrum of (|4.7p contains (2m + 1) singlet states, (2m — 1) of 
which correspond to bound states. The middle bound state and its energy are 

V'm,o(a;;6') = ^^^^j^™ , £:„(6')=mC0S^ , Hci^rafi{x--,0) = £'m{0)'4>rnfi{x]9) . (4.12) 

Note that the energy £m{d) = mcos| 7^ for 9 ^ tt, and that it varies monotonically from +m to —m when 9 
varies in the interval (0, 27r). For boundary values 9 = 0,2tt the bound state ipm,o{x;9) as well as all other bound 
states disappear, and the system has a spectrum of a free massive (1+1)D Dirac particle consisting of two symmetric 
scattering bands (—00, — m] and [m, 00). Thefore, for 6' 7^ {modir) the energy reflection symmetry of the real case is 
partially broken. 

Relations ()4.1ip and (|4.12p show that the middle bound states of the two BdG systems with parameter 9 shifted 
in 2ti satisfy the relation £m{9) = —£m{9 + 2tt). This equality together with (|4.1ip correspond to a kind of duality 
relation 

Hc{x;9) = -Hc{-x:9 + 2Tr) , (4.13) 

which is reflected in the structure of the spectrum shown on Fig. [31 cf. [ssj . 
Observe that 

(_ _ m(m-l)sin^ f _^ ^2 q \ 

dS^ cosh2(sinfx) I (A^A\ 

^ dx^ cosh2(sin|2^) + ™ / 

is diagonal, with Schrodinger-type Hamiltonians of Poschl- Teller form on the diagonal. In the general complex case 
of the BdG Hamiltonian, this local diagonalization of H^ is not possible. 

The non-trivial integral of motion which comes from the Lax pair relation of the nS+ hierarchy can be presented 
in the form 



_ ,2mt{-ircoslV'_Z-+\ {t>lZ)^' 

" \ Viz 

2mz(-l)" cos ^'D^jn+lt'-m+2--T^m-l X'-m+l^?-m+2 • • -I'm- l^^n 



(4.15) 



Since nonlocal operator TZa^ commutes with the Hamiltonian (|4.7p and anticommutes with S2rm 

[7^^T3,i^c] = o , {7^a3,52™}-o , (4.16) 

we identify the operators 

21=^2™ , Q2 = ^7^(73Ql (4.17) 
as the supercharges. They generate the N = 2 polynomial superalgebra of the form (|3.2ip . 

[Qa,i^c]=0 , {Qa,Qb}^25abPiUHc) , (4.18) 
in which the polynomial (|3.34p becomes 

2m- 1 

Pin.{Hc) = (He - £o)iHc - £2,n) J] (^c - £.f ■ (4.19) 
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FIG. 3: Spectrum of the BdG system (|4.7p for m = 2 and m = 3. The red lines correspond to the singlet states energies; the 
middle red line reveals a partial breaking of the energy reflection symmetry. The pairs of the same symbols on the vertical lines 
indicate singlet energies for dual systems with shifted in 2ti parameter 6. 



Let us discuss now the nature of supersymmetry of the assoeiated extended Schrodinger system described by the 
Hamiltonian (|4.14|) . Though (|4.14|) is a square of the BdG Hamiltonian (|4.7p . the latter can not be treated as a usual 
supercharge of linear Af — 2 supersymmetry with identified as a Z2-grading operator. The reason is that though 
the quadratic diagonal operator (|4.7p commutes with 173, the linear operator (j4.7p does not anticommute with the 
diagonal Pauli matrix. In the rest of this section we present three distinct but related possible interpretations of the 
supersymmetry associated with the second order matrix diagonal Hamiltonian 

Hc^Hl . (4.20) 

First, we note that the BdG Hamiltonian (|4.7p can be presented as 

Q 

Hc^ K-mcos-as , (4.21) 

where 

and 2?,„ is defined in (|4.8|) . Operator (|4.22p is the integral of motion for the system He- In correspondence with 
(|4.9p . it is the rescaled real BdG Hamiltonian Hk{x) defined in (|3.26p . H^{x) — sin |iJR(sin Unlike the BdG 
Hamiltonian (j4.2ip . it anticommutcs with (T3 and so, can be identified as one of two supercharges of the Af = 2 linear 
supersymmetry, Qi — H^. Defining the second supercharge in a usual way, Q2 ~ i<^3Qi, we arrive at the superalgebra 



{Qa,Qb} = 25,b (He- m^cos^f) , [Q„Hc]=0 



(4.23) 
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Since supercharges Qa are hcrmitian, from the first of these relations and relation (|4.20p we find that the energies of 
the finite-gap BdG system (|4.7p satisfy the inequality £^ > m? cos^ |. Since for 6* 7^ (morf 27r) system (|4.7p has an 
odd number of 2m +1 singlet states in its spectrum, the supersymmetric structure (|4.23p of the associated Schrodinger 
system implies that the reflection energy symmetry is broken for Q ^ ix (mod 27r) . 

As in the real case 9 = tt, relation (|4.20p means that the BdG nontrivial integral (|4.15p is also the integral of motion 
for the Schrodinger system Tic- Like the BdG Hamiltonian He, integral S2m does not have a definite Z2-parity with 
respect to the grading operator CT3, but analogously to (|4.21|) it also can be decomposed into bosonic (diagonal) and 
fermionic (antidiagonal) operators 

S2rn^J'''+B' . (4.24) 

JF^ is a rescaled nontrivial integral of a real case, T^{x) ~ (sin |)^™ S'B(sin |x), with Sti{x) denoting here the operator 
p.33p . Bosonic operator is 

S^ = 2m(-l)™cos^i^^™-i^n+ , (4.25) 

where 11+ = i(l + 0-3) is a projector. Both operators and are the integrals of motion for the Schrodinger system 
(|4.20p . Bosonic integral can be related to a rescaled operator (|3.1ip with index m changed for m — 1. One can 
find a superlagebra generated by the fermionic set of integrals Qa, Si = , Si = ia^Si, and by the bosonic integral 
Bg together with the Hamiltonian (|4.20p . It is a nonlinear extension of superalgebra (|4.23p . whose explicit form we 
do not display here. 

Instead of hermitian integrals Qa, a = 1,2, associated with decomposition ()4.21l) . consider now a pair of hcrmitian 
conjugate integrals 

which are certain linear combinations of the integrals Qi, Qi = ia^Qi and (T3. Together with the Schrodinger 
Hamiltonian Tic they generate the centrally extended Af = 2 superalgebra 

{Q,Q^^nc , [nc,Q]^[nc,Q^ = o , (4.27) 



{Q,Q} = {Q\Q^ = c , (4.28) 

in which the operator 

C = cos^ -I2 , (4.29) 

proportional to the 2x2 unit matrix I2, plays the role of the central charge. Relation Tic > C is implied by 
superalgebra (g^H), via relations (H^^f >0,H^ = i{Q - Qt). 

Central extension of A/" = 2 superalgebra of such a structure was discussed in [s^ within the framework of super- 
symmetric quantum mechanics. Since mutually conjugate matrix integrals (|4.26p do not anticommute with (T3, they 
can not be identified as supercharges if (T3 is treated as the Z2-grading operator T. In our case, however, we can 
identify one of the operators (|4.10p as F. They both commute with He, and identify the Schrodinger Hamiltonian 
He and the integral as the bosonic operators, while integrals (|4.26p are identified by any of them as the fermionic 
operators. Having in mind that [(73,(5] = —2iQ^, [o^3,Q^ = 2iQ, the operator cr^ can be treated then as a generator 
of i?-symmetry. 

Finally, let us discuss the tri-supersymmetric structure of the system (|4.20p . From the described properties of the 
finite-gap BdG system (|4.7p it follows that the involutive operators (|4.10p and their composition TZa^ are the trivial 
integrals of motion for the associated Schrodinger system (|4.20p . [Gi,Hc] — [R-cr3,Hc] — 0. Let us denote = He 
and = Sim, and introduce the operator 

= ]^{X\y'} . (4.30) 

Operators , and Z^ are nontrivial integrals for the Schrodinger system (|4.20p . which at = tt reduce to the 
integrals of the Schrodinger system corresponding to the real case. The choice of one of the three trivial integrals as 
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the Z2-grading operator classifies the nontrivial integrals as bosonic and fermionic operators in correspondence with 
relations 

[A'^7^a3]=0 , {y,7^(T3} = {Z^7^(73} = , (4.31) 

[y,7^T2^] = o , {A'^7^^2x} = {2^7^^2J = o , (4.32) 

[Z\azT2^] = Q , {X'',a3T2^} = {y'',a3T2^}^0 . (4.33) 

As we noted above, for 9 =^ tt the treatment of the operators 7?,T2^ and cr3r2ir requires a consideration of the pair 
of the corresponding dual systems. 

B. Complex periodic solution of the (2m — l)**" nS+ equation 

A periodic generalization of the hmction A(x) in ()4.ip that also satisfies the {2m — 1)*'' nS+ equation and includes 
as special case the examples (real and complex) discussed above is [the ra ~ 1 case was studied in 0]] 

A(x; 9) = exp M ^ C(^e/2) + * ns(*0/2)) + i 0%/2] , (4.34) 

where the parameter A sets the scale, A = A{9,v) = — 2i sc(i0/4) nd (i0/4). The functions a = (t(x; tJi, 0^2) and 
C = C(^i "-^ii ^2) are the Weierstrass sigma and zeta functions j53| . The real and imaginary half-periods are oji = K(i'), 
and a;3 = i K' = iK(l — v). The constant r/3 = C(iK') is purely imaginary. The parameter 9 plays the same role as 
in the previous non-periodic case, but now its period interval is [0, 8K'(!^)]. Using the properties of elliptic functions 
one can check that for the specific values of the parameter 9 = 2K' and 9 = 6K', (|4.34p reduces to the real periodic 
solution (|3.14p . In the infinite period limit ly ^ 1, the solution (|4.34p reduces to the complex non-periodic case (|4.ip . 
Meanwhile for 6* = 0, 6* = 4K' or = 8K', the solution reduces to a free particle case given by A = bie^^'^^ , where hi, 
1 = 1,2 are constants, that produces a simple shift in the energy 

To analyze the properties of function (|4.34p . we decompose it in terms of the amplitude and phase; A = Me*-*^, 



M{x) = |A| = mAyJV{Ax + iW) - V{i9/2) , (4.35) 

where 'P{x) is the elliptic Weierstrass function. Meanwhile M is a function of period T = 2'K/A in x, x changes in 
(f = 2K(-iC(i6'/2) + ms(i6'/2) - ^) over the period T. As a rcsuh ((TM)) changes as A(.t + 2K/A) e''^A(a;), , i. e. 
(|4.34p is a quasi-periodic in x function. Hence, in correspondence with l|3.ip . the BdG Hamiltonian (jl.ip satisfies the 
relation 

H{x + T) = e''<'''^/'^H{x)e-'^'-'fl'^ , (4.37) 

i. e. it is periodic up to a chiral transformation. This means that H{x + T) and H{x) systems are unitary equivalent, 
and the BdG Hamiltonian with quasi-periodic complex function ()4.34p describes the finite-gap periodic system. 

The spectrum of this system shown on Figure [¥] has the properties similar to those we described in the non-periodic 
case. The continuum bands (—00, Eq] and [E^m-i, 00), as well as any corresponding pair from 2m — 2 valence bands 
which do not include the middle band are symmetric with respect to the zero energy level. As in the non-periodic case, 
the spectrum reveals also a symmetry associated with a duality relation, which, up to a unitary chiral transformation, 
has the form of (|4.13p with 27r changed here for the 0-half-period equal to 4K'(i') [t^ . 

To reveal the Af = 2 nonlinear supersymmetric structure and the partial breaking of the energy reflection symmetry 
in the finite-gap BdG system given by (|4.34p . we identify first the Z2-grading operator. To this aim we note that the 
imaginary part of ()4.34p is an even function, while its real part is parity-odd. Applying the unitary transformation 
(|4.5p with parameter 6* = tt to the BdG Hamiltonian (|l.ip we obtain 

/ f(A-A*) *A + i(A* + A)\ / -ImA ^{£+RcA)\ 
''^-[^^^liA*+A) -i(A-A*) )'[^{£-RoA) ImA ) " ^'^'""^ 
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FIG. 4: Spectrum of the BdG system for the complex periodic solution (|4.34|) with m = 1 in dependence on 6. The pairs of the 
same symbols on the red lines corresponding to the singlet edge-states reflect a symmetry associated with the duality between 
the pairs of the systems with parameter 6 shifted in the half-period AK.' (y). 



The diagonal and anti-diagonal parts of the transformed BdG Hamiltonian have opposite parities, and similarly to 
the antiperiodic case, ()4.38p commutes with the nonlocal operator TZcr^, [T^-CTs, He] = 0. We can also find the analogs 
of the the operators (|4.3p which factorize the integral TZa^i and anticommute with (|4.38p . These are 

Qi = 'T3T4K' , Q2 = 7^^4K' , (4.39) 



{G^,Hc]^Q , i = l,2 . (4.40) 

Here we introduced the half-period ^-translation operator T^^i ■.9^9 + 4K'. 

As in the previous non-periodic case, from relation (j2.5p we can construct an operator 5*2™ which is a nontrivial 
integral of motion [He, S2m\ = 0. In contrast with the previous cases, here it is not easy to find a closed general form 
for this operator, and we just refer to Eq. (|2.5p . Conserved quantity S2m anti-commutcs with integral 7?.cr3, that 
allows us to identify the integrals Qi = S2m and Q2 ~ ilZcr^Qi as supercharges. They generate the nonlinear Af ^ 2 
supersymmetry, 

{Qa,Hc]=Q , {Qa,Qb} = 26abP4,n{Hc) , (4.41) 

where 

4m -1 

P4,n{Hc) = n (Hc-E,) , (4.42) 

i=l 

is the spectral polynomial, and Ei are the band-edge energies. In the limit ^ 1, (|4.42p reduces to the polynomial 
(|4.19p of the non-periodic case. 

The square of the BdG Hamiltonian (|4.38p takes here a non-diagonal form in contrast with the diagonal form of 
the associated Schrodinger Hamiltonian we had in the non-periodic case. Explicitly, we find that 

^ 2^/-^+Rc2A + (ReA)'+Im2A z(ImA)' A 

^ \ -i(ImA)' -^+Re2A-(ReA)' + Im2Ay 

includes purely imaginary off-diagonal terms. The second order Hamiltonian (j4.43p can be written in a more simple 
form if we define the matrix superpotential W, 

yy = ReAcr3-ImAcr2 , He = — r^ + W^ + W . (4.44) 



22 



Though the second order matrix Hamiltonian includes the ofF-diagonal terms, which can be interpreted as a coupling 
between the partner subsystems, the tri-supersymmetric structure still presents in this case as well. With the identi- 
fication of operators ~ He, ~ S2m and = ^{X,y}, as in the limit case (|4.30p when — > 1, we can treat 
them as nontrivial integrals of the second order system Tic- Their fermionic-bosonic nature and the concrete form of 
the corresponding superalgebra depend on which of the three trivial integrals 

7^^73 , 7^T4K' , or a3T4K' (4.45) 

is identified as the Z2-grading operator, cf. Eqs. (|4.3ip - (|4.33p . 

V. CONCLUSIONS 

The Bogoliubov/de Gennes Hamiltonian (|1.24p appears in a wide variety of contexts in theoretical and mathematical 
physics. Here we have shown first that the Ginzburg-Landau expansion of the associated grand canonical potential 
has a recursive structure that is precisely that of the AKNS integrable hierarchy. This explains why it is possible 
to solve the corresponding inhomogeneous gap equation by a simple ansatz for the Gorkov resolvent, as was done in 
[3]. We further have shown that this integrable AKNS hierarchy structure, and its finite gap solutions, provide the 
natural formalism for characterizing nonlinear quantum mechanical supersymmetry and tri-supersymmetry. We have 
illustrated these ideas with explicit solutions in terms of trigonometric and elliptic functions, for which the spectral 
and quantum mechanical supersymmetry properties can be computed in closed-form. Mathematically, there are open 
questions regarding the explicit form of the multi-dimensional theta function solutions to the AKNSat hierarchy 
equations for N > 2, which are associated with hyperelliptic [rather than elliptic] curves, and Riemann surfaces of 



J37I . [60l | . This is deeply related to the problem of the reduction of abelian integrals on Riemann 



genus g > 2 1131. 131.. .j. ^ ^ 

surfaces |3ll l37ll6Cll. 6ll |. An interesting case that is not yet well understood in such explicit terms is the complex AKNS 
solution with an even number of gaps. Physically, interesting open problems are: (i) the question of the existence of 
self-consistent solutions A(x) to the inhomogeneous gap equation for the massive NJL2 model, for which approximate 
solutions are known [gl], but no analytic solutions are known; (ii) the search for self-consistent crystalline solutions 
in higher dimensional models [i^, [gJ, [63|. It would also be interesting to apply these integrability and quantum 
mechanical sup ersymmetry structures to the generalized Gross-Ncveu and NJL models coming from string theories 

VI. APPENDIX 

In this Appendix we present the solution of the eigenvalue problem for the non-periodic BdG Hamiltonian (|4.7|) in 
terms of solution of the same problem for the real case 6 = tt p.26p . These two problems can be related in a simple 
way because of a very special nature of the BdG Hamiltonian corresponding to the twisted kink (j4.ip : in general case, 
in contrast with (j4.7[) . the square of (jl.l[) can not be transformed into a local diagonal form. 

Let us denote (|3.28p and (|3.31|) universally by -0+ and ip- as solutions with positive and negative energy respectively. 
As we noted, for zero energy eigenvalue, Eqs. (|3.28p and (|3.3ip give the same eigenstate p.30p . Due to the energy 
refiection symmetry of the real case, eigenfunctions of positive and negative energy sectors are related as ?/'_ = cr3'0+. 
We have 

H^{x)4)±{x) = ±sm^i+iij±{x) , '4)±{x) = i:±{sm ^x) , (6.1) 

where H^{x) = sin |77K(sin |x) corresponds to the notation introduced in (|4.2ip . and by we denoted the eigenvalue 
of the state Hk.'4'+ = £+ip+- Using relation (|4.2ip . we can rewrite (|6.ip as 

(_ffc + m cos Icts) sin |f +-0+ , {He + m cos ^a^) (73^+ = — sin |£+cr30+ . (6.2) 

From here we find the eigenvalues and eigenstates of the Hamiltonian He, 

Hc^P±ix;e)^±£+{0)i^±{x;0) , (6.3) 



£+{9) = xlsm^-£l+m^cos^- , ,p±{x; 9) = {a±{e) + a^H'+ix) , (6.4) 
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where 

a 

— 772 COS 

"+(^) = 1^7^^ ' a-(^?) = a+(0 + 27r) . (6.5) 

t+(t^) — sm 

Note that the middle bound state energy £m{0) = mcos |, for which £"+ = f+(7r) = in (|6.4p . is antiperiodic function 
of the parameter 9, £,n{S + Stt) = —£„i{9). For the rest of the spectrum, energy levels change periodically in 6 with 
the period 27r. that is reflected coherently in Fig. [31 
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